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ABSTRACT 

With  the  increasing  application  of  power  electronics  to  shipboard  electric  power 
systems,  there  is  a  need  for  studying  power  systems  in  an  object-oriented,  graphical 
environment  in  order  to  examine  the  interrelations  of  the  various  system  components. 

This  thesis  presents  a  method  of  analyzing  shipboard  electric  power  systems  in  a 
graphical  environment.  The  simulation  method  used  allows  for  real-time  observation  of 
the  various  power  system  measurements  of  interest 

The  machine  equations,  including  those  for  voltage,  torque,  and  rotor  angle  are 
developed  which  permit  analyzing  the  machine  in  any  chosen  reference  frame  depending 
on  the  study  being  undertaken.  Regulators  for  both  speed  and  field  voltage  are  developed 
and  the  stability  of  the  speed  regulator  is  examined.  A  means  of  modeling  a  load  for  a 
synchronous  generator  is  then  shown  which  accurately  interrelates  the  loading  of  the 
generator  and  the  frequency  and  voltage  output  of  the  machine.  This  load  is  then 
connected  to  the  synchronous  generator  and  two  different  scenarios  are  examined, 
including  a  fifty  percent  load  change  and  a  fault  study.  Graphs  showing  the  output  of 
these  studies  are  presented. 
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SYNCHRONOUS  MACHINE  MODEL  DEVELOPMENT 


The  primary  means  of  shipboard  electric  power  generation  is  with  synchronous 
generators  driven  by  steam  turbines.  The  electrical  and  electromechanical  behavior  of  the 
synchronous  generator  can  be  accurately  predicted  from  the  equations  for  a  three-phase 
salient-pole  synchronous  machine. 

The  rotor  of  the  synchronous  generator  is  equipped  with  a  field  winding  and  one  or 
more  damper  windings.  The  rotor  windings  usually  all  have  different  electrical 
characteristics.  Besides  the  differing  electrical  characteristics,  the  salient-pole  machine 
exhibits  a  magnetically  asymmetrical  geometry  due  to  the  varying  air-gap  path  for  the 
magnetic  flux.  This  magnetic  asymmetry  creates  time-varying  self-  and  mutual-inductance 
terms  in  the  voltage  equations  for  the  stator  variables.  By  transforming  the  equations  into 
a  reference  frame  which  rotates  at  the  same  angular  velocity  as  the  rotor,  the  magnetic  flux 
path  remains  constant,  thereby  eliminating  the  time-varying  inductance  terms. 

This  chapter  develops  the  voltage  and  electromagnetic  torque  equations  for  the 
synchronous  generator  in  terms  of  machine  variables  and  then  converts  these  equations  to 
the  rotor  reference  frame  via  Park's  transformation  [1].  The  developed  model  assumes  a 
magnetically  linear  system  thereby  neglecting  the  effects  of  saturation.  Saturation  could 
be  accounted  for  with  the  addition  of  another  simulation  segment  for  the  model,  however 
saturation  data  for  the  machines  studied  was  unavailable. 


A.     VOLTAGE  EQUATIONS  IN  MACHINE  VARIABLES 

Voltage  equations  for  the  synchronous  generator  will  be  developed  from  a  simple 

two-pole,  three-phase,  wye-connected,  salient-pole  synchronous  machine  shown  in  Figure 

1.   The  number  of  poles  of  the  particular  machine  can  be  easily  changed  as  this  is  one  of 

the  variables  of  the  resulting  torque  equations.  The  stator  windings  are  assumed  identical 

with  sinusoidally  distributed  windings,  displaced  120°,  with  Ns  equivalent  turns  and 

equivalent  resistance  rs  (subscript  s  denoting  stator  variables).  A  field  winding  and  three 

damper  windings  (fd,  kd,  kql,  and  kq2  respectively)  are  installed  on  the  rotor.     All 

windings  are  assumed  sinusoidal  with  applicable  N  equivalent  turns  and  resistance  r 

(subscript  denoting  particular  circuit).  The  quadrature  axis  (q  axis)  and  direct  axis  (d  axis) 

are  used  with  reference  to  the  magnetic  axis  of  the  rotor  circuits.   Voltages  of  the  rotor 
damper  windings  (vw,  v^p  and  v^2)    are  for  monitoring  purposes  only,  the  windings 

themselves  are  physically  shorted. 

Voltage  equations  in  machine  variables  [2]  may  be  expressed  in  matrix  form  as 


yaba~      rs1abcs^r  Probes 


(1) 

(2) 


where 
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Figure  1.  2-Pole,  3-Phase  Synchronous  Machine. 
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In  (1)  and  (2),  the  time  differential  operator  is  p.  In  (4)  and  (6),  the  subscript  /  denotes  a 
leakage  inductance  term.     The  subscripts  skql,  skq2,  sfd,  and  skd  denote  mutual 
inductances  between  stator  and  rotor  windings. 
The  magnetizing  inductances  are  defined  as 


Lm=-{LA-LB) 
Zw=|(LA  +  LB) 


(7) 
(8) 


where  the  term  incorporating  L4  refers  to  the  inductance  term  which  is  independent  of 
rotor  position  and  the  term  with  L#  is  the  inductance  term  which  is  position  dependent 


The  rotor  variables  are  then  referred  to  the  stator  windings  via  the  appropriate  turns 
ratio  giving  the  following  equations 
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where  the  primed  variables  (i.e.  L^.  )  refer  the  appropriate  rotor  variable  to  the  stator 
windings. 

The  voltage  equations  expressed  in  terms  of  machine  variables  referred  to  the  stator 
windings  are 


abcs 
1  qdr 


Ip(K)t    <+PK 


abcs 

i' 

lqdr 


(12) 


Note  that  the  voltage  equations  (12)  are  valid  for  positive  direction  of  stator  current 
assumed  out  of  the  stator  terminals.   If  motor  action  is  assumed,  the  sign  preceding  i^ 

must  be  changed. 


B.     TORQUE  EQUATION  IN  MACHINE  VARIABLES 

The  field  energy  stored  in  the  coupling  field  of  a  synchronous  machine  can  be 
expressed  as 


~    (1abcs)     (Ls    ~  LlsI)labcs  (labcs)     L^qdr 


Wf  =  -  (iabcsr  (L.  --  LbI)iabcs  -  uabcsy  ux 

(13) 

+  T(T)(iU)T(L;-L;i)i,qdr 


where 


U=diag[z^7     I^2    Ufd     l^\  (14) 

The  electromagnetic  torque  is  found  by  taking  the  partial  derivative  of  the  coenergy  with 
respect  to  the  rotor  angular  displacement  8r.  The  magnetic  system  is  assumed  linear, 
therefore  the  coenergy  equals  the  field  energy,  and  upon  taking  the  partial  derivative,  the 
electrical  torque  is 

Te  =  (|)|-|(i^)r^[L,-4l]i^  +  (UJr^[L;]^|  d5) 

where  P  is  the  number  of  poles  of  the  machine.  The  torque  is  positive  for  generator  action 
with  the  positive  direction  of  stator  current  assumed  out  of  the  stator  terminals. 
The  torque  and  rotor  speed  are  related  by 


Te=-J(j)P<»r  +  Tl  (16) 


where  J  is  the  inertia  of  the  machine  in  kilogram  meters2  (kg-m2)  or  Joule  seconds2  (J*s2). 
The  input  torque  Tj  is  positive  for  a  torque  applied  to  the  shaft  of  the  synchronous 
machine. 

C.  VOLTAGE  EQUATIONS  IN  THE  ROTOR  REFERENCE-FRAME 

The  machine  voltage  equations  (1)  and  (2)  are  seen  to  contain  flux-linkage  terms 
which  are  dependent  on  rotor  position  6r .  The  analysis  of  the  machine  behavior  is  greatly 

simplified  by  transforming  these  time-varying  inductances  to  constant  inductance  terms. 
This  can  be  realized  if  the  stator  variables  are  transformed  to  the  rotor  reference  frame. 
R.H.  Park  was  the  first  to  incorporate  a  change  of  variables  in  the  analysis  of  synchronous 
machines  [1]. 

A  change  of  variables  which  formulates  a  transformation  of  the  3-phase  variables  of 
the  stationary  circuit  elements  (thus  the  need  to  refer  the  rotor  variables  to  the  stator)  to 
the  rotor  reference  frame  may  be  expressed  as 
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The  term  K^  denotes  the  transformation  matrix  of  stator  variables  (the  subscript-^)  to  the 
rotor  reference  frame  (the  superscript-r).  It  can  be  shown  that  the  voltage  equations  (1) 
and  (2)  are  transformed  into 
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(22) 
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where 


V. 


ds 


ay  = 


-K 


qs 


0 


(24) 
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The  flux-linkages  can  be  expressed  in  the  rotor  reference  frame  by  transforming  (9)  to 
obtain 
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Expanding  the  voltage  equations  given  by  (22)  and  (23)  gives 
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(3D 
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and  the  expanded  flux  linkage  terms  become 
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Expressing  the  voltage  and  flux  linkage  equations  in  terms  of  reactances  rather  than 
inductances  gives  the  voltage  equations  in  matrix  format  (40),  where  the  expression  for 
the  applied  field  voltage  is  introduced 


<?,=v;^  (4D 

rfd 


and  the  leakage  and  magnetic  reactances  are  combined  as 

Xq  =  Xk  +  Xm  (42) 

Xd  =  Xb  +  X^  (43) 

Xkqx  =  XlkqX  +  X^  (44) 

%kq2  =  ^Ucq2  +  Xm  (45) 

X'fd  =  X;fd  +  Xw  (46) 

XL  =  X'm  +  X^  (47) 

The  zero  sequence  terms  have  been  omitted  since  balanced  conditions  are  assumed. 

The  model  being  developed  will  be  connected  to  a  non-infinite  bus.  This  implies  that 
the  voltage  and  frequency  do  indeed  vary  with  load,  based  on  the  droop  characteristics  of 
the  particular  regulator.  With  this  in  mind,  it  is  necessary  to  solve  for  currents  based  on  the 
bus  voltages,  therefore  the  currents  are  selected  as  independent  variables.  Rewriting  (40) 
with  the  current  and  current  derivative  terms  separated. 
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(48) 
The  state  equation  format  of  (48)  is  used  for  the  simulation,  so  the  current  derivatives 
must  be  solved  for.  As  an  example 


v  =  [A]i  +  [B]/?i 
pi  =  [B]-1{v-[A]i} 


(49) 


and,  with  this  format,  the  current  derivatives  are  solved  as  shown  in  (50) 
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Denoting  the  inverse  matrix  elements  as[B]      =[b]  = 


bn    &12 


'21 


'61 


the  desired  differential  current  equations  are  obtained 


til.  jr    _  m    ±md  yr   _  ±md  ZtT  n 

■qs  -W-~qs     ■    -s-qs     ■     — r   „     lds         wr     ,„       *«         "'r     „„       fta  J 

C06  C06  C06 


f)Z  =  b„[vr  +  rir  +  co  ^-Z  -  co 


'66 


(51) 


■^d   .•»•  ,*    ^mdif 


i  =  U4  +  rJL  +  ©r  -ft  5,  -  co 


co, 


co, 


*/B   "  fflr 


~^7rad  ,vr  -i 


CO, 


+  Wy£2   "  r^2^] 


(52) 


P^2  =  M£  +  r£  +  cor  ifi  -  ©r5*§  -  cor  i*C] 


CO, 


CO, 


CO, 


(53) 


p& 


CO*  CO,  CO. 


(54) 
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Pi'n  =  U<4  +  rfds  -  cor %l ?  +  cor  ^^  +  <or^i'^] 

C06  C06  C06 

+  W«3,  -  2^if]  (55) 

Pi*  =  Wii  +  rA  -  o^t  +  (Dr^SJ  +  cor  ^i- ] 

co6  co6  co6 

+  W«S,  "  2*&]  (56) 

Note  the  similar  form  of  the  q-  and  the  d-axis  integral  current  equations.    These  will 
greatly  simplify  the  SIMULINK  [3]  model. 

D.     TORQUE  EQUATION  IN  SUBSTITUTE  VARIABLES 

The  torque  equation  (15),  transformed  into  rotor  reference  frame  variables  is 

%  =  (f)(f  )^w(-&  +  i'i  +  %K  ~  W-&  +  Si  +  &WJ  (57) 

and  writing  in  terms  of  reactance 


^  =  — ^Zw=^-and4,=^-  (58) 

CO,  C0fc  CO,  ^         CO, 
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gives 


re  =  (l)A[iL(-/i +,- +%)rm  -^f-c-/;  +%  +i;r2)^]       (59) 

Z       Z        UK  COl 


't  wb 


The  rotor  speed  is  obtained  from  (16)  as 


w  =  (~27)(7;~:r/)  (60) 


E.     ROTOR  ANGLE 

The  rotor  angle  of  a  synchronous  machine  is  a  useful  parameter  to  monitor  the  torque 
production  of  the  machine.  It  relates  the  angular  displacement  of  the  rotor  and  stator 
magneto  motive  force  axes.  In  equation  form 

5  =  er  -  e„ 
=  JoWS)  -  <M$)W5  +  er(0)  -  ew(0)  (61) 

where  6  is  the  rotor  angle  (in  radians),  and  ^  is  a  dummy  variable  of  integration. 

Equation  (61)  can  be  implemented  by  solving  equation  (61)  in  terms  of  the  available 
parameters.  The  rotor  position,  9r  is  available  by  integrating  the  rotor  speed.  However, 
the  voltage  angle  0^  must  be  solved  for  as  follows 
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and 


vas  =  vqs  cos  6r  +  vA  sin  8, 


=  Jv^  +  vlcos[er  +  tan"1(^)] 

=  4vl  +  v*  cost  J  fflr  dt  +  0r  (0)  +  tan" l  (^)] 


v* 


(63) 


therefore,  from  (62) 


=  f©r<fr +  er(0)  +  tan-1(^-)  (64) 


er 


V 


giving 


5  =  9r-e„  =  tan"1^)  (65) 


where  a  positive  rotor  angle  is  assumed  for  generator  action. 
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H.         SPEED  AND  VOLTAGE  REGULATORS 


For  the  machine  model  to  accurately  simulate  a  shipboard  power  system  it  is 
necessary  to  include  controllers  for  both  the  field  voltage  and  the  prime  mover  torque. 
The  two  principal  controllers  of  a  turbine  generator  set  are  the  automatic  voltage  regulator 
and  the  speed  governor.  The  automatic  voltage  regulator  model  consists  of  voltage 
sensing  equipment,  comparators  and  amplifiers  controlling  a  synchronous  machine  which 
can  be  generating  or  motoring.  The  speed  governor  may  be  considered  to  have  similar 
equipment  but  in  addition  it  is  necessary  to  take  the  turbine  into  account 

A.     SPEED  REGULATOR 
1.     Model  Development 

The  block  diagram  of  a  steam  turbine  and  associated  governor  for  a  typical 
naval  ship's  service  turbine  generator  (SSTG)  is  shown  in  Figure  2  [4]. 

The  reference  speed  input  is  the  no-load  frequency  selected  by  the  operator  of 
the  electric  plant  The  gain  K  is  used  to  adjust  the  speed  droop  of  the  regulator.  zs  is  the 
time  constant  associated  with  the  servo  and  steam  valve  system  andzp  is  the  time  constant 
associated  with  the  conversion  of  steam  to  torque  (Tpm  ).  The  following  values  were  used 

for  this  study 

xs  0.1  sec 

%p  0.06  sec 

K  150 

From  this  block  diagram  (Figure  2),  a  SIMULINK  model  was  developed  and  is  shown  in 
Figure  3. 
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Figure  2.  Block  Diagram,  Steam  Turbine  &  Governor. 


19 


Figured  SIMULINK  Speed  Regulator. 
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As  can  be  seen,  the  speed  regulator  is  designed  to  compare  the  reference  speed  with  the 
speed  being  produced  by  the  synchronous  machine.  This  speed  difference  determines  the 
prime  mover  torque  being  applied  to  the  generator.  When  the  speed  difference  is  no 
longer  changing,  the  prime  mover  torque  is  equal  to  the  electrical  torque  developed  by  the 
generator  (load  dependent).  The  final  machine  speed  (and  therefore  the  machine  output 
frequency)  will  be  at  some  value  below  the  no-load  set-point  frequency.  This  can  be 
demonstrated  by  referring  to  the  'house'  diagram  shown  in  Figure  4. 

As  the  output  power  of  the  machine  is  increased,  the  operating  frequency  of  the 
machine  decreases  (based  on  the  slope  or  droop  of  the  particular  regulator)  the  operating 
frequency  can  be  restored  to  a  desired  frequency  (i.e.,  60  Hz)  by  the  operator  adjusting  the 
no-load  setpoint  of  the  regulator. 
2.     Regulator  Stability 

As  was  mentioned  in  the  previous  section,  the  speed  droop  characteristics  of  the 
regulator  are  determined  by  the  gain  K  of  the  regulator.  The  percent  speed  regulation  is 
defined  by 


co  ,-co„ 
%  Speed  Regulation  =  ~ L  ioo%  (62) 


As  the  gain  K  is  raised,  the  percent  speed  regulation  approaches  100%,  at  which 
point  the  regulator  and  therefore  the  regulated  plant  will  become  unstable.  This  stability 
point  was  determined  by  examining  the  root  locus  of  the  model  of  Figure  3.  The 
synchronous  machine  model  was  obtained  by  linearizing  the  machines  current  block 
around  its  stable  operating  point  for  rated  conditions  by  using  the  MATLAB  [5]  function 
linmod.  The  linearized  plant  model  returned  by  linmod  is  in  linearized  state  space  format 
which  is  then  changed  to  transfer  function  format  with  the  function  ss2tf.  The  transfer 
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Figure  4.  House  Diagram. 
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function  format  of  the  plant,  speed  block,  and  speed  regulator  is  then  reduced  to  a  single 
open  loop  to  study  the  roots  of  the  system  as  the  regulator  gain  K  is  varied.  This 
complete  process  is  demonstrated  in  Figures  5,  5a,  5b,  and  5c.  The  resulting  root  locus 
plot  is  shown  in  Figure  6.  The  operating  constraint  was  determined  for  the  first  root 
meeting  the  imaginary  axis  for  a  given  K.  The  limit  of  K  was  determined  to  be  153. 
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Figure  5.  Synchronous  Machine  Model. 
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Figure  5a.  Current  Block. 
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Figure  5b.  Linearized  Current  Block,  Speed  Regulator,  and  Speed  Block. 
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Figure  5c.  Transfer  Function  Development,  Speed  Regulator. 
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Figure  6.  Root  Locus  Plot-Synchronous  Machine  and  Speed  Regulator. 
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B.     VOLTAGE  REGULATOR 

The  voltage  regulator  or  excitation  system  used  in  this  study  is  shown  in  block 
diagram  form  in  Figure  7  [6].  The  parameters  are  defined  as  follows 

TA  and  KA  Regulator  amplifier  gain  and  time  constant 

T£  and  KE  Saturating  exciter  gain  and  time  constant 

Tn    '  XF2  » anc*  KF  Feedback  loop  gain  and  time  constants 

The  SIMULINK  model  developed  from  the  block  diagram  is  shown  in  Figure  8.  The 
regulator  compares  the  bus  voltage  with  the  specified  reference  voltage.  As  can  be  seen  in 
the  Simulink  model,  this  comparison  is  done  in  the  qdO  reference  frame. 
The  reference  voltage  is  computed  as  follows 


V   t  =i/V     +  V 

' reference       \     1  I 


—  2.-2 

-  1 


(63) 


This  equates  to  the  desired  line-to-neutral  voltage  (peak)  as  shown  below 


v      =-Kr\ 

V 

vd 

_2 
~3 

cos0r 
sin6r 

cos(6  -— ) 
r     3 

sin(9r-^) 

/r.     2tc. 
cos(0r  + — ) 
r      3 

2n 
sin(6r+— ) 
3 

V 

_vo_ 

1 

2 

1 
2 

1 
2 

yc_ 

(64) 


(65) 
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Figure  7.  Block  Diagram,  Voltage  Regulator. 


30 


^ 


Tut  tl 


eh      la 


Figure  8.  SDMULINK  Voltage  Regulator. 
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where 


va=vfacos(0J 

2% 
v6=vfacos(0CT-— -) 

2n 

Vc=VfaC0S(8CT+-y) 


(66) 


and 


and  solving  for  Vq  and  v^ 


er=£a)r®y+er(0) 
e„=J©.©y+e„(0) 


(67) 


v„  =— v. 


2ft  2ft  27C  2ft  " 

cose.cose^  +cos(0r  -— )cos(0cv  — — )+cos(0r  +— )cos(0cv  +— ) 


(68) 


Vd=TVi„ 


sinO.cos©^  +sin(0r  -  yOcos^  -y)+sin(0r  +y  ^os^  +  -y) 


(69) 


and  using  the  trigonometric  identities 


.      2iz.       .      2%.          ,      2tcx       ,      2tt.     3       ,        .  ,_«,. 

cosxcosy+cos(;t )cos(y )+cos(*+ — )cos(v+ — )=— cos(jt-y)  (70) 

~J  3  T~J  J  ±* 

.   ,      2tc.       .      2k.     .    .     2tcx       ,      2jc.    3  .   ,        .  ,_,,. 

sin  xcosy+sin(;c )cos(y )+sin(jc+ — )cos(y+ — )=— sin(jc-y)  (71) 

3  D  3  3  £ 
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gives 


v  =— v,  — cos(0  -0   ) 
q    3  h2         r      "  (72) 

=v/ncos[(cor-coe>+er(0)-eev(0)] 
v,=|v,|sin(er-ej  (?3) 

= vfaan[(o>,  -coe)f +0r(O)-0ev(O)] 


At  synchronous  speed  cor  =coe,  and  from  (61)  8=0r  -6W,  therefore 


and 


v9=vfacos(8)  (74) 

v(i=vi„sin(8)  (75) 


Reference    =    H  l™2  (5)   +   SUl*  (8)]    =    Vh  (76) 
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EDL       MACHINE  LOAD 


To  accurately  study  the  machine's  performance  it  is  necessary  to  be  able  to  load  the 
given  machine  in  a  realistic  manner.  The  inductive  nature  of  transmission  lines  and 
induction  motors  dictate  a  need  for  a  resistive-inductive  (RL)  loading  of  the  machine. 

This  section  develops  the  RL-load  used  for  the  synchronous  generator  and  then  the 
SIMULENK  implementation  of  this  loading. 

A.     PASSIVE  RL-LOAD 

The  load  used  for  this  study  was  comprised  of  a  passive  resistive-inductive  load. 
This  loading  is  typical  of  ac  motors,  transmission  lines,  electric  heating,  and  lighting.  The 
load  was  modeled  in  a&c-variables  in  order  to  eliminate  the  need  for  additional 
transformations  (abc  — >  qdO).  The  simulation  uses  voltage  as  the  input  and  currents  as  the 
outputs.  This  requires  the  phase  voltages  to  be  made  up  from  the  network  equations.  The 
loads  also  use  voltage  as  the  inputs  and  current  as  the  outputs.  Based  on  this  topology  it 
is  advantageous  to  assume  a  Y-connected  set  of  large  resistors  are  connected  across  the 
terminals  of  the  synchronous  generator  [71  and  the  phase  voltages  at  the  terminals  of  the 
generator  are  computed  from 

Va=  ~  WaG+iaL)  (85) 

Vr-%+iJ  (86) 

vc=-*OgHl)  (87) 

where  iaG,  i^,  and  icC  are  the  phase  currents  flowing  into  the  generator  from  the  load  and 

i^,  ibL,  and  icL  are  the  phase  currents  flowing  into  the  load  from  the  generator.   These 
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currents  are  computed  first  as  iq,  id,  and  i0  and  then  transformed  to  ia,  ib,  and  ie.   This 
is  represented  in  Figure  9. 

This  load  modeling  is  quite  accurate  as  the  resistance  value  R  is  made  large  (20-50 
pu).  As  the  resistance  is  made  large,  the  current  through  this  branch  becomes  negligible. 

B.     SIMULINK  LOAD  MODEL 

The  passive  RL-loa.d  was  modeled  in  SIMULINK  by  starting  with  the  voltage 
equations  (only  developed  for  phase  a,  others  similar)  from  Figure  10. 

va=ea+Lspia+rsia  (88) 

va=rL+LLpiL  (89) 

va=Wa~iL)  (90) 

The  differential  equations  can  be  implemented  as  shown  in  Figure  1 1. 

As  can  be  seen,  the  terminal  voltage  va  is  dependent  on  the  output  current  ia .  This 
voltage  dependency  is  the  difference  between  a  finite  and  an  infinite  bus  model  of  electric 
machines.  An  infinite  bus  analysis  assumes  no  dependency  between  va  and  ia . 

An  attempt  was  made  to  make  the  load  purely  resistive.  SIMULINK  is  unable  to 
deal  with  this  type  of  load  however  since  this  would  incorporate  direct  feed-through 
between  input  and  output  of  the  load.  This  forms  an  algebraic  loop  and  cannot  be 
analyzed.  The  inductance  of  the  transmission  line  overcomes  this  limitation.  The  purely 
resisitve  load,  connected  by  a  transmission  line  (with  inherent  inductance)  to  the  generator 
bus  adds  sufficient  inductance  to  warrant  modeling.  This  inductance  then  incorporates  a 
time-delay  between  input  and  output  due  to  the  integral  form  of  the  load  equations. 
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Figure  9.  Generator  Terminal  Voltage  Model. 
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Figure  10.  Passive  RL-Load  Model. 
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Figure  11.  Differential  Equation  Implementation  of  Voltage  Equations. 
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IV.       MACHINE  MODEL  PERFORMANCE 

The  synchronous  machines  performance  is  examined  in  this  section.  The  machines 
parameters  are  given  and  the  machines  response  to  both  a  load  change  and  a  fault 
condition  are  studied. 

A.     MACHINE  PARAMETERS 

The  machine  being  simulated  is  a  4-pole,  60  Hz,  440  V  (line-to-line),  500  kW,  0.85 
power  factor  synchronous  generator.  The  per  unit  parameters  on  the  machine  base  (Z  ^^e 
=  0.1291  Q)  are  as  follows 

stator  resistance,  rs  =  0.00031 

rotor  field  resistance,  /yy  =  0.00226 

q-axis  damper  winding  resistance,  r^,  =  0.0224 

d-axis  damper  winding  resistance,  r^/  =  0.00673 

stator  leakage  reactance,  Xfc  -  0.0632 

d-axis  magnetizing  reactance,  X^j  =  1.4648 

q-axis  magnetizing  reactance,  X^  =  0.7668 

field  leakage  reactance,  X/#/  =  0.1545 

q-axis  damper  winding  leakage  reactance,  X/^  =  0.0771 

d-axis  damper  winding  leakage  reactance,  X/&/  =  0.3676 

rotor  inertia,  J  =  14.526 
These  machine  parameters  are  taken  from  the  Machinery  Systems  Laboratory  (MSL)  test 
facility  for  a  typical  shipboard  synchronous  generator  [8].  The  MSL  test  facility  is  within 
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the  Machinery  Systems  Engineering  Branch  at  David  Taylor  Naval  Ship  Research  and 
Development  Center. 

B.     SIMULATION:  LOAD  CHANGE 

The  machine  was  initially  loaded  at  an  approximate  power  factor  of  0.8,  with  the  load 
calculated  as  follows 

Machine  data:      440-V    1.5  MVA .  From  Figure  12 


11  3      254 

/  =  1.969  e-Ja*~%m)  kA 


Z  = 


254  e 


jo 


=  l29eiXM°mQ 


(91) 
(92) 
(93) 


1.969*f;36-87  kA 
rL  =  103. 2  mQ.  XL  =  77.42  mQ  (94) 

With  the  machine  loaded  as  above,  the  expected  steady-state  power  delivered  to  the 
load  will  be  approximately  1.2  MW.  This  does  not  however  account  for  the  losses  in  the 
stator  windings.    The  machine  then  underwent  a  50  %  step  load  change  by  doubling  the 
values  of  r .  and  Xf  .  This  was  done  in  SEMULINK  by  the  load-change  block  shown  in 
Figure  13.  The  plots  of  the  response  are  shown  in  Figures  14  through  21.  The  data 
obtained  from  the  simulation  was  as  follows: 

FULLY  LOADED 

MODEL 

1.1754 

1.1792 

59.99 

0.7992 


Load  power,  MW 
Machine  power,  MW 
Frequency,  Hz 
Power  Factor 


50%  LOADED 
MODEL 

0.5186 

0.5151 

64.32 

0.7832 
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Figure  12.  Initial  Load  Parameter  Determination  Circuit 


41 


Figure  13.  SEVfULENK  Load  Change  Block. 
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Figure  14.  50  Percent  Load  Decrease,  vq. 
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vd,  50%  load  decrease 


200- 


>100 


0 


0 


0.05 


0.1 
time  (sec) 


0.15 


0.2 


200- 


>100 


0 


■ 

1        '  -   1 

i 

■ 

i 

0 


2  3 

time  (sec) 


Figure  15.  50  Percent  Load  Decrease,  vd. 
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va,  50%  load  decrease 
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Figure  16.  50  Percent  Load  Decrease,  va. 
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Figure  17.  50  Percent  Load  Decrease,  field  voltage. 
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Figure  18.  50  Percent  Load  Decrease,  ia. 
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te,  50%  load  decrease 
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Figure  19.  50  Percent  Load  Decrease,  te. 
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del,  50%  load  decrease 
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Figure  20.  50  Percent  Load  Decrease,  rotor  angle  (del). 
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wr,  50%  load  decrease 
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Figure  21.  50  Percent  Load  Decrease,  rotor  speed  (a^). 


50 


C.     SIMULATION:  THREE-PHASE  FAULT 

The  stability  of  a  synchronous  machine  throughout  a  power  system  following  a  fault 
is  of  major  concern.  A  3-phase  fault  or  short-circuit  rarely  occurs  and  a  3-phase  fault  at 
the  machine  terminals  is  even  more  uncommon;  however,  it  is  still  instructive  to 
demonstrate  the  ability  to  insert  a  fault  for  a  chosen  time  interval  and  study  the  machines 
dynamic  performance. 

The  machine  was  initially  loaded  at  a  power  factor  of  0.8  as  in  part  B.  A  3-phase 
fault  was  inserted  and  was  left  in  place  for  7  cycles.  The  fault  impedance  (representative 
of  the  bus  work)  is  0.001291  +  j  0.1291  Q.  If  the  fault  had  been  allowed  to  remain  on  the 
system  slightly  longer  the  machine  would  have  become  unstable  in  terms  of  not  having 
returned  to  synchronous  speed  after  removal  of  the  fault  or  the  machine  would  have 
slipped  poles  prior  to  returning  to  synchronous  speed.  It  was  determined  that  a  seven- 
cycle  fault  is  representative  of  a  typical  breaker  reclosure  timing  cycle  which  the  machine 
could  have  been  subjected  to. 

The  three-phase  fault  was  performed  in  SEMULINK  by  the  fault  block  shown  in 
Figure  22.  The  fault  implemented  by  using  a  step  function  block  which  stepped  in  the 
fault  impedance  at  a  particular  time  and  then  another  step  function  block  stepped  in  the 
normal  load  at  the  end  of  the  desired  fault  time.  The  plots  of  the  response  are  shown  in 
Figures  23  through  30. 
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Figure  22.  SEVfULINK  Three-Phase  Fault  Block. 
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Figure  23.  7-Cycle,  3-Phase  Fault,  vq. 
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Figure  24.  7-Cyde,  3-Phase  Fault,  vd. 
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Figure  25.  7-CycIe,  3-Phase  Fault,  va. 
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Figure  26.  7-Cycle,  3-Phase  Fault,  field  voltage. 
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Figure  27.  7-Cycle,  3-Phase  Fault,  ia. 
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Figure  28.  7-Cycle,  3-Phase  Fault,  te. 
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Figure  29.  7-Cyde,  3-Phase  Fault,  rotor  angle  (del). 
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wr,  7-cycle  3-phase  fault 
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Figure  30. 7-Cyde,  3-Phase  Fault,  rotor  speed  Ceo,.). 
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V.         CONCLUSIONS 


A.  ADVANTAGES 

The  advantages  to  studying  naval  power  system  components  in  a  graphical 
environment  are  many.  To  be  able  to  see,  in  real  time,  any  line  or  load  parameter  (i.e. 
current  or  voltage)  makes  the  study  of  system  component  interaction  very  straightforward. 

The  concept  of  "closing  the  loop"  to  make  the  input  and  output  of  the  machine 
dependent  on  one  another,  as  opposed  to  the  infinite  bus  concept,  is  easily  implemented 
with  the  graphical  tools  available.  Navy  power  systems  definitely  require  this  non-infinite 
bus  modeling  due  to  their  capacity  difference  as  compared  with  commercial  power  grids. 

B.  ANALYSIS  OF  SIMULINK 

Matlab's  SIMULINK  is  an  excellent  tool  for  the  type  of  studies  conducted  in  this 
thesis.  Model  building  and  integration  into  the  "power  system"  was  easily  done. 

One  of  the  difficulties  with  SIMULINK  included  the  algebraic  loop  problem 
necessitating  an  integrator  in  the  path  to  provide  a  time  delay  element  In  this  particular 
study  this  was  not  a  problem  due  to  the  resulting  differential  equations,  however,  a 
possible  purely  resistive  load  would  be  troublesome. 

On  the  issue  of  speed  of  the  simulations,  naturally  one  is  always  looking  for  a  faster 
simulation.  I  began  my  research  on  a  486-33MHz  machine,  then  graduated  to  a  SPARC-2 
machine  and  realized  a  speed  difference  of  about  two  times.  I  then  moved  to  a  SPARC-10 
and  achieved  a  speed  increase  of  approximately  four  times  over  the  SPARC-2.  Typical 
times  to  run  a  simulation  of  approximately  5-seconds  was  about  4-hours  on  the  SPARC- 
10.  In  comparing  this  to  a  non-graphical  simulation  tool  such  as  Advanced  Continuous 
Simulation  Language  (ACSL)  [9],  a  comparable  model  runs  approximately  in  real  time  (60 


61 


seconds  of  run  time  data  in  approximately  60  seconds).  The  disadvantage  of  ACSL  is,  of 
course,  the  inability  of  the  user  to  observe  system  output  in  real  time  and  the  non-modular 
modeling  of  systems. 

C.         FUTURE  WORK 

The  modeling  of  additional  load  components  including  induction  machines,  another 
synchronous  generator  sharing  the  load,  parallel  RL-loading,  and  adjustable  speed  drives 
could  all  be  performed  in  SIMULINK.  As  the  number  of  integrators  in  the  complete 
model  increases,  however,  the  run  time  will  also  increase.  In  addition,  the  existing  RL- 
load  can  be  examined  under  such  conditions  as  critical  clearing  time,  first  swing  transient 
stability  limit,  and  various  fault  conditions  such  as  single  line  to  ground  faults,  line-to-line 
faults,  etc. 

In  addition,  voltage  regulator  stability  can  be  studied  to  find  gain  limits  for  this 
particular  regulator.  Different  regulator  designs  could  also  be  examined. 

Mathworks  has  also  advertised  a  software  package  which  will  take  existing 
SIMULINK  models  and  compile  them;  this  will  supposedly  increase  simulation  speed 
considerably.  This  product  sounds  very  encouraging. 


62 


LIST  OF  REFERENCES 

[1]        Park,  R.H.  "Two-Reaction  Theory  of  Synchronous  Machines-Generalized 
Method  of  Analysis-Part  L"  AIEE  Trans.,  Vol.  48,  July  1929,  pp.  7 16-727. 

[2]        Krause,  Paul  C,  Analysis  of  Electric  Machinery,  McGraw-Hill,  Inc.,  1987. 

[3]        SIMULINK,  Dynamic  System  Simulation  Software,  The  Math  Works,  Inc.,  Natick, 
Mass.,  February  1993 

[4]        Krause,  Paul  C.  and  Wasynczuk,  O.,  "Development  of  Computer  Models  to 
Analyze  Shipboard  Electrical  Power  Systems,"  Final  Report,  Contract  No. 
N61533-89-D-0037,  School  of  Electrical  Engineering,  Purdue  University,  West 
Lafayette,  IN,  October  1990 

[5]        MATLAB,  High-Performance  Numeric  Computation  and  Visualization  Software, 
The  Mathworks,  Inc.,  Natick,  MA.,  February  1993 

[6]        Arrillaga,  J.,  Arnold,  C.P.,  Computer  Analysis  of  Power  Systems,  John  Wiley  & 
Sons  Ltd.,  1990. 

[7]        Krause,  P.C.,  Final  Report  on  "Modeling  of  Shipboard  Electric  Power  Distribution 
System,"  for  David  Taylor  Research  Center,  SBIR  Phase  I  Contract  No.  N00167- 
88-C-0031,  P.C.  Krause  and  Associates,  Inc.,  West  Lafayette,  IN,  July  1988 

[8]        Mayer,  C.B.,  "Digital  Computer  Simulation  Program  for  Turbine  Emulator  Drive 
and  Load  System,"  Report  of  Contract  N00167-85-C-0131.  May  13, 1987. 

[9]        Advance  Continuous  Simulation  Language  (ACSL),  Mitchell  and  Gauthier 
Associates  (MGA)  Inc.,  Concord  MA.,  1993 


63 


BIBLIOGRAPHY 


Anderson  P.M.  and  Fouad  A.A.,  Power  System  Control  and  Stability,  Iowa  State 
University  Press,  1977. 

Glover  J.D.  and  Sarma  M,  Power  System  Analysis  and  Design,  PWS-Kent  Publishing 
Company,  1989. 

Sarma,  M.,  Synchronous  Machines,  Their  Theory,  Stability,  and  Excitation  Systems, 
Gordon  and  Breach  Science,  1979. 

Stevenson,  W.D.,  Elements  of  Power  System  Analysis,  McGraw-Hill,  Inc.,  1982 


64 


INITIAL  DISTRIBUTION  LIST 

1 .  Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  Virginia     22304-6145 

2.  Library,  Code  52 

Naval  Postgraduate  School 
Monterey,  California     93943-5002 

3.  Chairman,  Code  EC 

Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 
Monterey,  California     93943-5121 

4.  Professor  Robert  W.  Ashton,  Code  EC/Ah 
Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 

Monterey,  California     93943-5 1 2 1 

5.  Professor  Harold  A.  Titus,  Code  EC/Ts 
Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 

Monterey,  California     93943-5121 

6.  Professor  Stephen  Williams 
Department  of  Engineering 
Baylor  University 

PO  Box  97356  " 

Waco,  Texas  79798-7356 

7.  LT  Daniel  Kline 
3708  Yates  Ave.  No. 
Crystal,  Minnesota   55422 

8.  Jack  Stowe 

6670  Vernon  Ave.  S.  #110 
Edina,  Minnesota   55436 


65 


DUDLEY  KNOX  udkaRY 
NAVAL  POSTGRADUATE  SCHOOi 
MONTEREY  Cm  ^3943-5101 


GAYLORD  S 


